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Abstract 

The Sagdeev potential technique has been employed to study the dust ion acoustic solitary waves 
and double layers in an unmagnetized collisionless dusty plasma consisting of negatively charged 
static dust grains, adiabatic warm ions, and isothermally distributed electrons and positrons. A 
computational scheme has been developed to draw the qualitatively different compositional param¬ 
eter spaces or solution spaces showing the nature of existence of different solitary structures with 
respect to any parameter of the present plasma system. The qualitatively distinct solution spaces 
give the overall scenario regarding the existence of different solitary structures. The present system 
supports both positive and negative potential double layers. The negative potential double layer 
always restricts the occurrence of negative potential solitary waves, i.e., any sequence of negative 
potential solitary waves having monotonically increasing amplitude converges to a negative poten¬ 
tial double layer. However, there exists a parameter regime for which the positive potential double 
layer is unable to restrict the occurrence of positive potential solitary waves. As a result, in this 
region of the parameter space, there exist solitary waves after the formation of positive potential 
double layer, i.e., positive potential supersolitons have been observed. But the amplitudes of these 
supersolitons are bounded. A general theory for the existence of bounded supersolitons has been 
discussed analytically by imposing the restrictions on the Mach number. For any small value of 
positron concentration, there is no effect of very hot positrons on the dust ion acoustic solitary 
structures. The qualitatively different solution spaces are capable of producing new results for the 
formation of solitary structures. 
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I. INTRODUCTION 


Theoretical and experimental stndies of the nonlinear dynamics of Ion Aconstic (lA) 
waves have received a great deal of attention for several decades. Using a mechanical analogy, 


Sagdeev established that nonlinear lA waves can exist in the form of periodic or solitary 
waves. On the other hand, in the same year, Washimi and Taninti nsed rednctive 
pertnrbation method to investigate the small amplitnde ion aconstic solitary waves in a 
Dseqnently, the nonlinear theory of lA waves was developed by a nnmber 
Experimental resn 


cold plasma. Sn 

n” 

of anthors p 


15 


ts of Ikezi et al. 


16L Nakamnra and Tsnkabayashi 


17l |. Nakamnra and Tsnkabayashi Nakamnra et al. 19|, Nishida and Nagasawa 20| . 


Nakamnra 


2l| and Cooney et al. j22| conhrmed the existence of lA waves. 


There has been considerable interest in stndying ion aconstic solitons and donble layers 
in electron-positron-ion (e-p-i) plasmas as snch plasmas are fonnd in supernovas, pulsar 
environments, cluster explosions, active galactic nuclei etc. Popel et al. |23| considered the 
nonlinear propagation of lA waves in cold plasma consisting of cold ions, and isothermally 
distributed electrons and positrons. They have found the existence of compressive solitons 
only. In this paper, they reported that the presence of the positron can result in reduction 
of the lA soliton amplitudes. Subsequently, the properties of the lA solitary structures in 
different electron-positron-ion (e-p-i) plasmas have been investigated by a number of authors 


24 


48|. 


However, in astrophysical environments highly (negative/positive) charged micronsize 
impurities or dust particulates are observed in addition to earlier mentioned e-p-i plasma. 
The presence of dust grains having large masses introduces several new aspects in the prop¬ 


erties of the nonlinear waves and coherent structures 


49 


59|. Depending on different time 
scales, there can exist two or more acoustic waves in a typical dusty plasma. Dust Acoustic 
(DA) and Dust Ion Acoustic (DIA) waves are two such acoustic waves. Shukla and Silin [50 1 
were the hrst to show that due to the quasi-neutrality condition rieo + naoZd = Uio and the 
strong inequality rieo nn) (where neo,njo, and n^o are, respectively, the number density 
of electrons, ions, and dust particles, and is the number of electrons residing on the 
dust grain surface), a dusty plasma (with negatively charged static dust grains) supports 
low-frequency Dust Ion Acoustic (DIA) waves with phase velocity much smaller (larger) 
than electron (ion) thermal velocity. In the case of a long wavelength limit, the dispersion 
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relation of DIA wave is similar to that of lA wave for a plasma with Ueo = Uio and Ti <^Te, 
where Ti{Te) is the average ion (electron) temperature. Due to the usual dusty plasma ap¬ 
proximations (rzeo <C Uio and Tj ~ Te), a dusty plasma cannot support the usual lA waves, 
but a dusty plasma can support the DIA waves of Shukla and Silin [5^. Thus, DIA waves 
are basically lA waves modihed W the presence of heavy dust particulates. The theoretical 


prediction of Shukla and Silin 
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50| was supported by a number of laboratory experiments 


6l| . The nonlinear theory of DIA waves in different du^ plasma systems has been 


investigated by Bharuthram and Shukl a 16211 , Nakamura et a l. 

651 , Shukla and Mamun jdol . Verheest et al. jddi . Sayed and Mamun [6 


and Shukla 


63| , Luo et al. [6^ , Mamun 


Aline]ad 68|. Baluku et al. 69|] . Das et al. 0 - 

It is well known that a fully ionized universe contains electrons, positrons, ions and 
micron-sized charged dust grains {m, 72|. The presence of e-p-i-d plasma has been detected 
in active galactic nuclei, pulsar magnetospheres, inter-stellar clouds, supernova environments 
as well as in laboratory experiments of cluster explosions by intense laser beams 


73 


78|. 


Therefore the nonlinear dynamics of the propagation of waves in this plasma system 
emerged as an interesting held to explore. For the hrst time, Ghosh and Bharuthram 


has 


79| 


investigated the nonlinear propagation of small but hnite amplitude ion acoustic solitons and 
double layers in a collisionless unmagnetized e-p-i-d plasma consisting of cold ions, negatively 
charged static dust particulates and Boltzmann distributed electrons and positrons. They 
have used the reductive perturbation method to derive the Korteweg-de Vries (KdV) equa¬ 
tion and the modihed KdV equations for the ion acoustic waves in the e-p-i-d plasma system. 
Using KdV equation, they proved that this e-p-i-d plasma system supports both compressive 
and rarefactive solitons if the parameters of the system satisfy certain conditions. Using the 
modihed KdV equation, they have derived the necessary conditions for the existence of weak 


double layers. Employing the reductive perturbation technique, El-Tantawy et al. 


80| inves¬ 


tigated the ion acoustic solitary structures in a collisionless unmagnetized four-component 
e-p-i-d plasma consisting of warm ions, superthermal electrons and positrons, and negatively 
charged static dust impurities. Using Bernoulli’s pseudo-potential method, Dubinov et al. 


8l| elaborated the nonlinear theory of dust ion acoustic waves in a collisionless unmag¬ 


netized four-component e-p-i-d plasma consisting of warm ions, isothermal electrons and 
positrons, and negatively charged static dust impurities. They have found the existence of 
four diherent types of physically realizable nonlinear ion acoustic waves: a subsonic periodic 
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wave, a large amplitude subsonic periodic wave (superlinear wave), supersonic rarefactive 
and compressive solitary waves. 


In the above mentioned papers of nonlinear DIA or lA waves in e-p-i-d plasmas, the 
authors have not considered the following interesting points. From the investigations of 
DIA waves in a collisionless unmagnetized three-component e-i-d plasma consisting of warm 
ions, isothermal or nonthermal electrons and negatively charged static dust impurities, we 


have seen from the article of Das et al. 


70| that the system supports the Negative Potential 


Solitary Waves (NPSWs) for all M > Me, where Me is the lower bound of the Mach number 
M for the existence of DIA solitary structures, i.e., DIA solitary structures start to exist for 


all M > Me- From this simple conclusion of Das et al. 


70|, we have the following questions 


regarding the existence of different DIA solitary structures when an amount of positrons is 
injected in the system: Can we restrict the occurrence of NPSWs after injecting positron 
in the system? Does the system support Negative Potential Double Layer (NPDL)? Does 
the NPSW disappear from the system if the concentration of positron increases? Does the 
system support Positive Potential Double Layer (PPDL)? Does the system support the co¬ 
existence of double layers of both polarities? Does the system support supersolitons? (i.e.. 
Is there any soliton after the formation of double layer? ) Is the amplitude of the super- 
soliton bounded? What is the condition on the Mach number for the existence of bounded 
supersolitons (supersolitons of bounded amplitude)? or. Can we impose any restriction on 
the Mach number to get bounded supersolitons? Is there any connection between the occur¬ 
rence of suppersolitons and the occurrence of NPSWs and / or Positive Potential Solitary 
Waves (PPSWs)? Does there exist a critical value of positron concentration for which the 
system supports only positive potential solitary structures (PPSWs and PPDLs)? Does 
there exist a critical value of positron concentration for which the system supports only 
Positive Potential Solitary Waves (PPSWs)? What happens if very hot positrons are in¬ 
jected in the system? or Is there any qualitative difference in the DIA solitary structures if 
very hot positrons are injected in the system? To answer all these questions, we reconsider 
the same problem of Das et al. 70| by making the following changes: instead of taking 
three component e-i-d plasma, a collisionless unmagnetized four component e-p-i-d plasma 
has been considered in which electrons and positrons both are isothermally distributed, he., 
in the present investigation, we reconsider the problem of Ghosh and Bharuthram [79j] in 
the following directions: (i) the equation of pressure for ion fluid is taken into account to 
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include the effect of ion temperature, (ii) instead of considering the reductive perturbation 
method, Sagdeev [l| pseudo-potential technique has been employed to investigate arbitrary 
amplitude DIA solitary waves, double layers and supersolitons in unmagnetized collisionless 
e-p-i-d plasma, (hi) the investigations have been made with the help of the qualitatively 
different compositional parameter spaces instead of considering particular values of the pa¬ 
rameters involved in the system, (iv) a general theory has been given on the existence of 
bounded supersolitons which describes the condition on the Mach number for the occurrence 
of bounded supersolitons, (v) a thorough investigation on existence of DIA solitary struc¬ 
tures has been presented through the compositional parameter spaces, giving the special 
emphasis on the existence of bounded supersolitons. 

Four basic parameters of the present e-p-i-d plasma system are p, p, ape and ate which 
are, respectively, the ratio of unperturbed number density of positrons to that of the total 
unperturbed number density of positive charge, the ratio of unperturbed number density of 
electrons to that of the total unperturbed number density of negative charge, the ratio of the 
average temperature of positrons to that of electrons and the ratio of average temperature 
of ions to that of electrons. The studies of the DIA solitary structures of the present e-p-i-d 
plasma system have been made over the entire physically admissible values of p, p, ape and 


The present paper is organized as follows: In IJTTl the basic equations are given. The 
derivation and the mechanical analogy of the energy integral are given in l lllll Conditions 
for the existence of solitary wave and double layer solutions are also given in this section. 
A general theory for the formation of bounded supersolitons has been presented in §IVI by 
imposing the restrictions on the Mach number. The lower bound and the upper bound of the 
Mach number for the existence of different solitary structures have been determined in this 
section. With the help of the analytical theory discussed in l lIVl a computational scheme 
has been developed in ^to draw the qualitatively different compositional parameter spaces 
with respect to any parameter of the system. In l|Vl DIA solitary structures of the present 
plasma system have been thoroughly presented with the help of the qualitatively distinct 
compositional parameter spaces, giving special emphasis on the bounded positive potential 
supersolitons. Finally, a brief summary and discussions have been given in 1 1 VII 
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II. BASIC EQUATIONS 


The following are the governing equations describing the non-linear behaviour of dust 
ion acoustic waves propagating along x-axis in collisionless unmagnetized dusty plasma con¬ 
sisting of adiabatic warm ions, negatively charged immobile dust grains, and isothermally 
distributed electrons and positrons: 

3ti ■ 3 

^ + ( 1 ) 


(dui dui\ dpi d(j) 

= 0’ 

\ ot OX J ox Ox 


( 2 ) 


dpi 

dt 


, dpi dui 

Ox Ox 


0 , 


(3) 


(9^0 

dx"^ 


= -ine^rii - Ue + Up - Z^Ud). 


(4) 


Here n,, rze, np, Ud, Ui, pi, 0, x and t are, respectively, ion number density, electron number 
density, positron number density, dust particle number density, ion fluid velocity, ion fluid 
pressure, electrostatic potential, spatial variable and time, 7 (= 3) is the adiabatic index, rrii 
is the mass of ion fluid, Zd is the number of negative unit charges residing on the dust grain 
surface and e is the charge of an electron. 

The above four equations are supplemented by 


Up 


^eO 


= exp 



Up 

0 


5 - — c-X p 

crpe^_ 


(5) 


and the charge neutrality condition 


^^0 “1“ ^eO “1“ 


( 6 ) 


where rZeo, rijO; and n^o are, respectively, the unperturbed number densities of electron, 
ion, positron and dust particulate, <h and ape are given by 
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(7) 


Here is the Boltzmann constant, Te and Tp are the average temperatures of electrons 
and positrons respectively. 
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III. ENERGY INTEGRAL 


The linear dispersion relation of the DIA wave for the present dusty plasma system can 
be written as 


= Cd\ I — ^ ^ 9 ) CD = CsMs 


( 8 ) 


k V 1 + k^xi 

where uj and k are respectively the wave frequency and wave number of the plane wave 
perturbation, and 


c. = M. =. hcTi, + 

m V p + iiape 


(9) 
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ClQ — UiQ + UpQ — UeO + Z^n^f). 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


Now for long-wave length plane wave perturbation, i.e., for /c —x 0, from linear dispersion 
relation ([8]), we have. 


lim — = cd and lim — = cd 
k^o k k^o dk 


(14) 


and consequently the dispersion relation ([8]) shows that the linearized velocity of the DIA 
wave in the present plasma system is cd with Ad as the Debye length. 

To study the arbitrary amplitude time independent DIA solitary waves and double layers, 
we make all the dependent variables depend only on a single variable ^ = x — Ut where U 
is independent of x and t. Thus, in the wave frame moving with a constant velocity U the 
equations ([I])-(lll) can be put in the following form 


4|(c/-„,)i9 ^ = 0, 

[ riio 


(15) 


^ ( TTo, _L -L <^2 f I ^ ^2 ^*0 d , Pis 


(16) 
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( 17 ) 


- l^ie 1 - ^ ri^ p Up Tlj 

Af, 1-p l-pUeO l-pUpO Tlio 


(18) 


Here P = riioKsTi, Uio is the unperturbed ion number density and Tj is the average tem¬ 
perature of ions. 

Using the boundary conditions, 


( —, Ui, 0, ^ (1,1, 0, 0, 0) as 1^1 ^ cx) 

nio r ai; 


( 19 ) 


and solving fflSj) . ffTbjl . and ffT7)l . we get a quadratic equation for n'f, and the solution of the 
hnal equation of Uj can be put in the following form: 


where 


a = jv, = (t^/g.)V2 

Tlio !^ I ! ^!ij ^ 

Y $ $ + Y 
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+ VScTie) 
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<h 


21^“ vScTie) 


( 20 ) 


( 21 ) 


Now integrating flTSll with respect to cj) and using the boundary conditions flT^ . we get 
the following equation known as energy integral with W{(j)) as the Sagdeev potential or 
pseudo-potential: 

Y |)% ir (^)=0 ( 22 ) 

where 




IV. + - p 

1 — p 1 — p 
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fUe = exp 



l,lUp 


1 




(25) 
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The energy integral (l2^ can be regarded as the equation of energy of a particle of unit mass 
moving along a straight line whose position is (j) at time ^ with velocity The hrst term 
of the energy integral fl22ll can be regarded as the kinetic energy of the particle whereas the 
second term of the energy integral (j2^ can be regarded as the potential energy of the same 
particle at that instant. As the kinetic energy is a non-negative quantity, W{(j)) < 0 for the 
entire motion of the particle. Differentiating the energy integral fl22|) with respect to 0, we 
get the following equation 


d?(j) 


+ 1D'(0) 


O,1T'(0) 


dW 

dcf) 


(27) 


The above equation shows that the particle of unit mass is under the action of the force 
and this force is attracting one, i.e., directed towards the point 0 = 0, i.e., 0 = 0 
is the centre of the force if hT'(0) > 0 for 0 > 0 or W'{(j)) < 0 for 0 < 0. Now, if 
1T(0) = IT'(O) = 0, then both the velocity and the force acting on the particle at 0 = 0 are 
simultaneously equal to zero and consequently, the particle is in equilibrium at 0 = 0. Let 
0 = 0 be an equilibrium position of the particle associated with the energy integral fl2^ . 
The nature of the equilibrium point 0 = 0 depends on the sign of 1T"(0). In fact, we have 
the following cases depending on the sign of Wifi) : Wfi > 0 and Wfi) < 0. The case 
W'fi) > 0 implies that 0 = 0 is a position of stable equilibrium but this case is not possible 
because for this case the radius of curvature of the curve hT(0) at 0 = 0 is positive and 
consequently the curve W (0) is concave with respect to the 0 - axis in a neighbourhood of 
0 = 0 which again implies that W (0) > 0 for any 0 (except the point 0 = 0) lying within a 
very small neighbourhood of 0 = 0. But W (0) > 0 is not possible for the particle which has 
been associated with the energy integral fl2^ as we have seen earlier that W (0) < 0 for the 
entire motion of the particle. Next we consider the case W'fi) < 0. For this case, 0 = 0 is 
a position of unstable equilibrium of the particle and W (0) is convex with respect to the 0 - 
axis in a neighbourhood of 0 = 0 which again implies that W (0) < 0 for any 0 (except the 
point 0 = 0) lying within a very small neighbourhood of 0 = 0. So, if the particle placed 
at 0 = 0 be slightly displaced towards the positive (negative) direction of 0 - axis, it moves 
away from its position of unstable equilibrium and it continues its motion until its velocity 
is equal to zero, i.e., until 0 takes the value 0m such that Wfim) = 0 for 0m > 0 (0m < 0). 
Now, if hF'(0m) > 0 for 0m > 0 iW'fim) < 0 for 0m < 0), then although the velocity of the 
particle at 0 = 0m vanishes, the force acting on the particle at 0 = 0m is directed towards 
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the point 0 = 0 and consequently, the particle will come back again at 0 = 0. Therefore, 
if hh(0) = IT'(O) = 0, 1T"(0) < 0, = 0 and > 0 for 0^ > 0 (1T'(0™) < 0 

for (prn < 0) then the energy integral (122]) can be regarded as the equation of energy of the 
oscillatory motion of a particle of unit mass, i.e., we have an oscillation of the particle within 
0 < 0 < 0m (0m < 0 < 0). As tliis oscillatiou takes place in a wave frame moving with a 
velocity U, this oscillation propagates with velocity U and in this wave frame this oscillation 
is known as positive (negative) potential soliton. From the mechanical analogy of the energy 
integral (122|) as mentioned above, Sagdeev ^ established that for the existence of a Positive 
Potential Solitary Wave (Negative Potential Solitary Wave) solution of the energy integral 
fl22p . the following three conditions must be simultaneously satisfied. 

Sa :: 0 = 0 is the position of unstable equilibrium of the particle, i.e., W{0) = W'{0) = 0 
and W”{0) < 0. 

Sb :: kF(0m) = 0, kF'(0m) > 0 for some (pm > 0 (kF'(0m) < 0 for some (pm < 0). 
This condition is responsible for the oscillation of the particle within the interval 
min{0, 0m} < 0 < max{0, 0m}. 

Sc :: W{(p) <0 for all 0 < 0 < 0m (kF(0) < 0 for all (pm < (p < 0). This condition is 
necessary to dehne the energy integral fl22|) within the interval min{0, (pm} < 0 < 

max{O,0m}. 

On the other hand, if W'{(pm) = 0 along with the condition W{(pm) = 0 then the velocity 
and the force acting on the particle at 0 = 0m are simultaneously equal to zero and conse¬ 
quently the particle will not be reflected back again at 0 = 0. In this case instead of soliton 
solution Energy Integral fl22|) gives shock-like solution which is known as double layer solu¬ 
tion. If (pm > 0 the double layer is known as positive potential double layer (PPDL) whereas 
if 0m < 0 the double layer is known as negative potential double layer (NPDL). Therefore, 
for the existence of a Positive Potential Double Layer (Negative Potential Double Layer) 
solution of the energy integral (1221) . the following three conditions must be simultaneously 
satished. 

Da :: 0 = 0 is the position of unstable equilibrium of the particle, i.e., W{0) = W'{0) = 0 
and W”{0) < 0. 
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Db :: iy(0m) = 0, iy'(0m) = 0, W'\(l)m) < 0 for some 0m > 0 (0m < 0). This condition 
actually states that the particle cannot be reflected back again at 0 = 0. 


Dc :: 1T(0) < 0 for all 0 < 0 < 0m (1T(0) < 0 for all 0m < 0 < 0). This condition 
is necessary to dehne the energy integral fl22l) within the interval min{0, 0m} < 0 < 
max{O,0m}. 


Here, we have not discussed the case for the existence of solitary structures (solitons and 
double layers) when 1T"(0) = 0. Following Das et al. 82| one can discuss this case for this 
particular problem. 

Therefore, the necessary conditions for the existence of solitary structures of the energy 
integral ([22]) are W{0) = 0, hF'(O) = 0, and hF"(0) < 0. It can be easily checked that 
hF(0) = 0, Wifi) = 0, and the condition fF"(0) < 0 gives U > cd, and consequently the 
greatest lower bound felb) of the Mach number (M = U/cd, which has been explicitly 
discussed by Dubinov |83|) for the existence of solitary structures of the energy integral fl2^ 
is Me = 1, i.e., the solitary structures start to exist for U > M > = 1). So we 

have introduced the following dimensionless quantities: x = x/Xd, ^ = ^/Xd, t = t/^Xo/ci)), 
Wi = Ui{XD/ cd)IXd = UijcD-, 0 = 0/*h, ^ = Pi/P, nj = ng/rio. Now we note the following 
fact: ^ = x — Ut = XjXx — U^t = Xd[x — = Xd[x — Mt\ ^ = x — Mt ^ = x — Ml, 

i.e., here spatial coordinate is normalized hy Xd and the time is normalized by ■^. Then with 
respect to these dimensionless quantities, the energy integral can be simplihed as follows, 
where we drop overline on both independent and dependent variables: 


1 

2 



+ V{ct>) = 0, 


(28) 


where 


H(0) = (Mf - 3am) 


V) + 


P 


1 — p 






1 — P 


Ve - \-^Vd 
1 — P 


Vi = M^M/ + am - Ni{M^M/ + Sa^ - 20 - 2amiVf), 

^ ^ __ MMs\/2 _ 

^io iV^M — 0 + — 0 ) 


(29) 


(30) 
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(mMs -f x/ScTje'j , 

(31) 

d'M = -| 

i^MMs - y/3ai^ , 

(32) 

Ve = e^-l,Vp = l- = 0. 

(33) 


The equations (1^ and (1^ are dynamically equivalent and consequently, the qualitative 
behaviour of solitary structures determined by these two equations are same. In fact, one 
can take any equation equivalent to the energy integral to discuss the qualitative be¬ 
haviour of the solitary structures, but in order to avoid quantitative inaccuracies, the general 
prescription is to normalize length by A/j and time by Ad/c_d, where A/? is the actual Debye 


length and co is the actual acoustic speed of the plasma system [Dubinov 83|]1. One can 
obtain cd by considering the linear dispersion relation as prescribed by Dubinov |83|. So, we 
can take the energy integral (j28ll to investigate the solitary structures of the present plasma 
system. 

From the expression of Ni as given by (l30|) . we see that Ni is well - defined if and only 
if both — 0 and — 0 are real and non-negative. Now the conditions — 4> ^ 0 
and Tm — 0 > 0 hold simultaneously if and only if 0 < Min{^M, As < d’jw, the 
equation (jHni) gives theoretically valid expression of Ni if and only if 0 < Tm- Therefore, 0 
is restricted by the inequality: 0 < Tm = {MMg — \/3^)^/2. In fact, we have used this 
inequality, viz., 0 < Tm to find the upper bound of the Mach number for the existence of 
PPSWs. 


IV. ANALYTICAL THEORY FOR THE FORMATION OF SUPERSOLITONS 

In this section, our main aim is to hnd the Mach number for the existence of supersolitons. 
Before going to investigate the existence of supersolitons analytically, hrst of all, we note 
the following facts: 

For the existence of solitary structures of the energy integral fl25D . we have D(0) = 0, 
D'(0) = 0, and 1A"(0) < 0. It can be easily checked that D(0) = 0, D'(0) = 0, and the 
condition l/"(0) < 0 gives M > Me = 1. Therefore, the solitary structures start to exist for 
M > Me = 1, i.e., M = Me = 1 is the lower bound of the Mach number M for the existence 
of solitary structures. 
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If there exist solitons after the formation of double layer of same polarity for increasing 
values of the Mach number then the solitons are known as supersolitons. In other words, 
suppose we have a double layer solution at M = of some polarity and for dehniteness, 
let us assume that we have a positive potential double layer solution at M = M^l. If 
there exists a Mach number Mp{> Mpp) such that we have a PPSW at M = Mp{> 
Mpi), then the PPSW at M = Mp{> Mpp) is called a positive potential supersoliton. 
Similarly, one can dehne negative potential supersolitons. Existence of supersolitons for 
some particular values of the parameters involved in different plasma system have been 
recent^reported by many authors 3, 84~99|. In the above investigations, except Das 
et al. [701], the authors have studied the existence of supersolitons for some particular values 


of the parameters involved in different plasma systems but Das et al. 


70| have studied 


the unbounded negative potential supersolitons through the compositional parameter space. 
They have clearly pointed out the existence region of the unbounded negative potential 
supersolitons but instead of ‘supersolitons’ they have used the term ‘dias type solitons’. 


From the physical interpretation of the energy integral for the existence of the double 
layer of any polarity, we have found that the existence of the double layer of any polarity 
always implies that there must exist a sequence of solitary waves of same polarity having 
monotonically increasing amplitude converging to the double layer solution, i.e., the ampli¬ 
tude of the double layer solution acts as an exact upper bound or least upper bound {lub) 
of the amplitudes of at least one sequence of solitary waves of same polarity. Therefore, if 
the double layer solution exists then this double layer solution is a limiting structure of a 
sequence of solitary waves of same polarity. On the other hand, any sequence of solitary 
waves ends with a double layer of same polarity if it exists. So, double layer solution plays 
an important role to restrict the occurrence of at least one sequence of solitary waves of 
same polarity. More specihcally, if M = Mpppp (M = M^ppi) corresponds to the positive 
(negative) potential double layer then there must exist positive (negative) potential solitary 
waves for any M restricted by the inequality Me < M < Mpppp {Me < M < Mnpdl) and 
the amplitude of the solitary wave increases with increasing M and these solitary waves end 
with a double layer solution of same polarity when M assumes the value M = Mpppi for 
positive potential double layer (M = M^ppi for negative potential double layer). Therefore, 
double layer solution plays an important role to restrict the occurrence of solitary waves of 
same polarity. Of course, our aim is to hnd a Mach number Mp > Mpppi (Mjv > M^ppp) 
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such that there exists a PPSW (NPSW) at M = Mp (M = Mn). Now, if the plasma 
system does not support any PPDL (NPDL) for any given set of values of the parameters 
of the system, then there is no question of existence of positive (negative) potential super- 
solitons. So, it is important to investigate whether the system supports any double layer 
solution. Following Das et al. [^, we shall analytically investigate the existence of double 
layer solution of the energy integral fl28l) . 

For the existence of a double layer solution of the energy integral fl28|) . we must have a 
non-zero 0 (0 7 ^ 0 ) such that the following conditions are simultaneously satished: 

1/(0)=O,D'(0) = O,1/"(0)<O (34) 


Using equations fl2^ - fl33ll . the hrst equation, the second equation, and the third inequality 
of (134)) can be written, respectively, as 


V(^) = (1 - Ni)AfAf, + - Ni(3a^ -2^- 2ai,Nf ) -S = 0, (35) 


V"(^) = iV, - ^ = 0 , 


where 


d 


dS\ 




5 = 




-K-f 


1- fi 


V.- 


pCTpe 


K, 


(l_p)- ■ (i-p)'“ (i-p) ^ 

Eliminating Ni from fl3^ and (l36|) . we get 


= 


S-^^e+% 


3aie - 20 - 2aie ( g 


2 n 




(36) 

(37) 


(38) 


(39) 


It can be easily checked that 0 = 0 if and only if dS/d(j) = 1, and consequently for non-zero 
0 , fl3^ can be written as 


a-Pa-P, = H^), 


(40) 


where 


h(0) = 




CTi, 


+ % [3cri, 


20 - 2cri, 


\ dtp ) J 


(1 


d4>) 


(41) 
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Using dUD, from fl36|) and (l37|) we, respectively, get 


= 0 . 


( 42 ) 



(43) 


where 


h(0) 


dS 


dS 


(44) 


\j9+{4>) -0 + -0 df 


^±(0) = ± \/3^)^ 


(45) 


Now, the double layer solution of the energy integral fl28l) having amplitude \4>di\ exists at 
M = Mdi, where Mdi is given by the following equations 



(46) 


V(<t>di) = 0 , 


(47) 


with 


h(M - > 0 , 


(48) 


9-{4>dl) — 4>dl > 0 , 


(49) 


drj 

*-*di 


(50) 


To derive condition (09]), we have used the following restriction on 0 : 0 < Tm = \{MMs — 
■v/ScTje)^. If the inequalities given by f08|) . (0^ and fl50|) hold simultaneously then one can 
get a PPDL or NPDL at 0 = (pdi according to whether (pdi > 0 or 0^; < 0. 

So, if the system supports any double layer solution, we can easily find the Mach number 
corresponding to that double layer solution by using the equations (0^ - (071) and the 
conditions (051) - (1501) . Now, if there exist solitary waves after the formation of double 
layer then our aim is to investigate whether the occurrence of PPSWs (NPSWs) after the 
formation of PPDL (NPDL) are bounded or not. The amplitudes of PPSWs or NPSWs are 
also monotonically increasing for increasing values of Mach number after the formation of 
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double layer. Now we want to investigate whether there exist any Mach number for which 
the solitary wave takes the largest amplitude after the formation of double layer. If such 
Mach number exists then this Mach number restricts the occurrence of PPSWs (NPSWs) 
after the formation of PPDL (NPDL) and consequently the existence region of PPSWs 
(NPSWs) becomes bounded. 

Secondly, if the system does not support either PPDL or NPDL or both then our aim is 
to investigate whether the occurrence of PPSWs or NPSWs are bounded or not. So, in this 
case our task is to hnd the upper bound of the Mach number (if exists) for the occurrence 
of PPSWs or NPSWs. In the next two subsections, we shall analytically investigate the 
upper bounds of the Mach number for the existence of PPSWs (including positive potential 
supersolitons) and NPSWs (including negative potential supersolitons) of the energy integral 


A. Upper bound of the Mach number for the existence of PPSWs (including 
positive potential supersolitons) 


To hnd an upper bound of the Mach number for the existence of any type of positive 
potential solitary structures, following Das et al. [70|, we consider the existence of a PPSW 
of amplitude 0m > 0. Then, following conditions are simultaneously satished. 


1/(0) = 0, U'(0) = 0 and U"(0) < 0 (51) 

l^(0m) = 0 and U'(0m) > 0 (52) 

U(0) <0 for all 0 < 0 < 0m (53) 


The conditions as given in (|5Til are simultaneously satished if M > Me = 1. We assume 
that this condition holds good. We also assume that the conditions as given in fl5^ are also 
true. Consider the condition as given in ([53]). 

We have seen in the previous section that D(0) is real only when 0 < Tm and conse¬ 
quently, we must have 0m < , otherwise is not a real quantity. Therefore, we 

can rewrite the inequality as given in fl53D as 


U(0) < 0 for all 0 < 0 < 0m < ^m (54) 

But the inequality flSTl) can dehne a large amplitude PPSW of amplitude which is in 
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conformity with (l5TD if 


V{^m) = 0 and V'{^m) > 0. (55) 

Again, let M^ax be the maximum value of M up to which positive potential solitary wave 
solution can exist. As 'hjvr increases with M, Therefore, the inequality 

V{(j)) <0 for all 0 < 0 < (0™ < Tm (56) 

dehnes the largest amplitude PPSW if 

V = 0 and V'(^M^ax) > 0- (57) 

Therefore, for the existence of PPSWs, the Mach number M is restricted by the following 
inequality: Me < M < Mmax, where M^ax is the largest positive root of equation V (Tm) = 0 
subject to the condition Vi^^M) > 0 for all M < Mmax- 

So, Mmax can restrict the existence of PPSWs. But if a PPDL exists at M = Mppdl, 
then this PPDL can restrict the existence of at least one sequence of PPSWs because the 
existence of a PPDL implies the existence of at least one sequence of PPSWs which converges 
to the PPDL solution. Therefore, with respect to the existence of Mpppi and Mmax, we 
have the following three cases. 

Case-1 : li Mpppi exists but Mmax does not exist for hxed values of the parameters 
involved in the system then obviously Mpppi is the upper bound of M for the existence of 
positive potential solitary waves, i.e., one can get a PPSW for all M such that Me < M < 
MpppL and for M > Mpppp there does not exist any positive potential solitary structure 
whereas at M = Mpppi one can get a PPDL solution. 

Case-2 : If Mmax exists but Mpppi does not exist for hxed values of the parameters 
involved in the system then Mmax is the upper bound of M for the existence of positive 
potential solitary waves, i.e., one can get a PPSW for all M such that Me < M < Mmax 
and for M > Mmax there does not exist any positive potential solitary structure. 

Case-3 : If both Mmax and Mpppi exist for hxed values of the parameters involved in 
the system then the case for which Me < Mmax < Mpppi is not possible for the following 
reasons: The existence of PPSW is restricted by the inequality Me < M < Mpppi only when 
there does not exist Mmax- In other words, existence of Mmax implies that the mach number 
M for existence of any type of positive potential solitary structure is always restricted by the 
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inequality Me < M < Mmax- In this case, we have seen that M^ax exists hnitely but Mmax 
is not the upper bound of the mach number M for the existence of PPSWs, which is not 
possible because there exist PPSWs for all M lying within the interval M^ax < M < Mpp£,L. 
So, if Mmax exists hnitely, then Mmax niust be the upper bound of the mach number M for 
the existence of any type of PPSWs, i.e., there does not exist any positive potential solitary 
structures for M > Mmax- 

Therefore, if both Mmax and Mpppp exist for hxed values of the parameters involved in 
the system then we must have M^ < Mpppp < Mmax- For this case we have the following 
conclusions: For this case, we can split the entire range of the Mach number M into two 
disjoint subintervals, viz.. Me < M < Mpppp and Mpppp < M < Mmax- Now, for 
Me < M < MppDL, we get a sequence of PPSWs converging to the positive potential double 
layer solution at M = Mpppp- In other words, the positive potential double layer solution 
at M = Mpppp can restrict the occurrence of all PPSWs for all M lying within the interval 
Me < M < MpppL whereas for Mpppp < M < Mmax, we get positive potential solitons 
after the formation of positive potential double layer at M = Mpppp, i-e., we get dias type 
solitons (according to Das et al. [^) or supersolitons (according to Dubinov and Kolotkov 


87l|). which has been recently reported by many authors for some particular values of the 
parameters involved in the system. For Mpppp < M < Mmax, amplitude of the supersoliton 
increases with increasing values of M but the amplitude of the supersolitons are bounded 
and the maximum amplitude of the supersoliton can be obtained at M = Mmax- So, we 
have two types of PPSWs: the hrst type is restricted by Me < M < Mpppr. whereas the 


second type is restricted by Mj 


PPDL 


< M < Mmax- In fact. Das et al. 


70| clearly stated 


the following conditions for existence of supersolitons or dias type solitons. Supersolitons 
can exist if there exist two types of solitary waves of same polarity separated by a double 
layer of same polarity and in this case there exists a jump type discontinuity between the 
amplitudes of the solitary waves separated by a double layer. In their paper, they have also 
reported that if there exists two types of NPSWs (PPSWs) separated by a NPDL (PPDL) 
then there is a hnite jump between the amplitudes of two types of NPSWs (PPSWs) only 
when dV/dM < 0 for all M > 0 and all 0 < 0(0 > 0). 


On the other hand, if there exists a PPDL solution 0 = 0 ppz)l(> 0) of the energy 
integral ( 128 |) for the Mach number M = Mpppp, then 0 = 4>ppdl is the smallest posi¬ 
tive double root of the equation D(0) = 0(= D(Mpppip,0) = 0) such that 0 < (j)ppDL < 
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^MppDL i-e-, V{MppdlAppdl) = 0 and V'{Mppdl,(1)ppdl) = 0 along with the condition 
0 < (t)ppDL < "^Mppop- If 0 = 4>ppdl{> 0) is the only root (donble root) of the eqnation 
l/(0)(= V{Mpppp,(j)) = 0) then this PPDL solntion is the nltimate solntion of the energy 
integral fl25]) and in this case, it is not possible to get any other PPSWs for M > Mpppp, 
i.e., for the occnrrence of any PPSW, the Mach nnmber M is restricted by the ineqnality 
Me < M < MppDp. Bnt if there exists a simple positive root 0 di(> 0) of V{Mppdl, 0) = 0 
for the nnknown 0 jnst after the occurrence of the double layer solution at 0 = (j)ppDL{> 0) 
for the Mach number M = Mpppp, i.e., if V{MpppL,(j)p)i) = 0 along with the condi¬ 
tions V{Mpppp, (ppppi) = 0, V'{Mppp)i, (ppppi) = 0, 0 < (j)ppDL < 4 >d\ < ^Mppop and 
V{Mpppp,(f)) < 0 for all (ppppL < 0 < 0z)i, then there exists a PPSW solution of the en¬ 
ergy integral fl25]) for at least one value of M > Mpppp, and consequently, we get a PPSW 
solution of the energy integral fl28|) after the formation of PPDL at M = Mpppp, i.e., there 
exists positive potential supersoliton. Although the amplitudes of supersolitons are mono- 
tonically increasing with increasing values of M for M > Mpppp, but the amplitudes of 
the supersolitons are restricted by since 0 < 0 < 'Lm- Therefore, one can hnd a Mach 
number M for which the amplitude of the supersolitons attains its maximum value and 
from the previous discussion, we have seen that the amplitude of the supersolitons attains 
its maximum value when M = Mmax{> Mpppi). 

B. Upper bounds of the Mach number for the existence of NPSWs (including 
negative potential supersolitons) 

We have seen earlier that D(0) is real if 0 < Tm , where Tm is strictly positive. For 
NPSWs or NPDLs, we have 

U(0) < 0 for all (pm < 0 < 0 (58) 

along with the conditions stated in Section UTTI for the existence of NPSWs or NPDLs. As Tm 
is strictly positive and for NPSWs or NPDLs 0 < 0, the condition 0 < is automatically 
satished and consequently for these two cases (NPSWs and NPDLs) D(0) is well dehned 
for all 0 < 0 without imposing any extra condition. Since there is no such restriction on 
0, we cannot use the same dehnition as in the case of PPSWs to hnd the upper bound of 
Mach numbers for the existence of NPSWs. For the case of NPSWs, to hnd an upper limit 
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or upper bound of M up to which NPSW can exist, we shall hrst hud a value M^pdl of 
M for which energy integral fl28|l gives a NPDL solution at M = M^pdl with amplitude 
0 = 4>npdl- Now, if at M = 0 = (Pnpdl is the only root (double root) of equation 

1/(0) = 1/(M, 0) = 0, i.e. V{Mj^p£)p, (pj^ppp) = 0 and V'^Mj^ppp, (j)j^p£)p) = 0, then the 
NPDL solution is the ultimate solution of the energy integral fl28p . and in this case no 
NPSW solution can be obtained for M > M^ppp, i-e. for the occurrence of NPSWs the 
Mach number M is restricted by the inequality Me < M < M^ppp. 

On the other hand, if there exists an inaccessible simple negative root 0 _d 2(< 0) of 
V{Mpfppp, 0) = 0 for the unknown 0 such that \(j)D2\ > \(I>npdl\ , i-e., iiV{MNppp, (f)p2) = 0 
along with the conditions V{Mnpdl,4>npdl) = 0, V'{Mnpdl ,4>npdl) = 0, |0d2| > 
\4>npdl\, and V^Mjpppp, (j)) < 0 for all 0 d2 < 0 < (I^npdl, then there exists a NPSW 
solution of the energy integral fl^ for at least one value of M > Mmpdl- Therefore, the 
double layer solution is unable to restrict the occurrence of NPSWs or there exist NPSWs for 
all M > MjsipdL) i-e., in this case, there exist negative potential unbounded supersolitons. 
For example, one can consider the DIA solitary structures in the plasma system considered 


by Das et al. 


70|. 


On the basis of the analytical theory, which have been considered in the present section, 
we have developed a computational scheme to draw the solution space or the compositional 
parameter space showing the nature of existence of different solitary structures of the present 
system. More specihcally, a solution space or compositional parameter space is a hgure where 
we plot the curves M = M^, M = Mmax, M = Mpppp and M = Mpfppp with respect to 
any parameter of the system and consequently, in the compositional parameter space, if 
we can hnd a region where the existence region of PPSWs (NPSWs) is separated by the 
curve M = Mpppp (M = M^ppp), then we can claim that the system supports positive 
(negative) potential supersolitons. In particular, for the case of PPSWs, if we can hnd a 
region where Mpppp < Mmax then we can claim that the system supports positive potential 
supersolitons for all Mach number M lying within the interval Mpppp < M < Mmax- So, if 
we can draw the solution space then it is simple to hnd the existence region of supersolitons 
and we can avoid the ‘trial and error’ method to investigate the existence of supersolitons 
for particular values of the parameters involved in the system. In the next section, we have 
investigated diherent solitary structures associated with the diherent solutions of the energy 
integral fl28|) with the help of the qualitatively diherent compositional parameter spaces. 
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V. DIFFERENT SOLUTION SPACES OF THE ENERGY INTEGRAL 


Four basic parameters of the present e-p-i-d plasma system are p, /i, ape and die, which 
are, respectively, the ratio of unperturbed number density of positrons to that of the total 
unperturbed number density of positive charge, the ratio of unperturbed number density 
of electrons to that of the total unperturbed number density of negative charge, the ratio 
of the average temperature of positrons to that of electrons and the ratio of the average 
temperature of ions to that of electrons. As the parameter aie assumes a constant value for 
dusty plasma system, we will discuss the existence regions of the different solitary structures 
with respect to p for different values of p and ape and with respect to p for different values of 
fi and ape- A solution space or compositional parameter space is a hgure where we plot the 
curves M = M = Mmax, M = Mppdl and M = Mnpdl with respect to any parameter 
of the system. In the present paper the solution spaces have been drawn with respect to the 
parameter p or p . To interpret the compositional parameter space we have made a general 
description as follows: 

1. Solitary structures start to exist just above the lower curve M = Me = 1. 

2. At each point on the curve M = Mpppp (M = M^ppp), one can get a PPDL (NPDL) 
solution. 

3. In absence of Mpppi M^ax (Mpppi) is the upper bound of M for the ex¬ 

istence of PPSWs i.e there does not exist any PPSW if M > M^ax {M > Mpppp). 
Although, it is important to note that there exist PPSWs along the curve M = Mmax 
but there does not exist any PPSW along the curve M = Mpppi. 

4. If both MppDL and Mmax exist hnitely, then ma.x{Mmax, Mpppp] is the upper bound 
of M for the existence of PPSWs. The case has already been discussed in subsection 
IIV Al of section [IVl 

5. If we pick a. p or p and go vertically upwards, then all intermediate values of M 
bounded by the curves M = Me and M = Mmax or Mpppp or niax{Mmax-i Mpppi} 
would give PPSWs. 

6. Similarly, all intermediate values of M bounded by the curves M = Me and M = 
Mnpdl would give NPSWs. 
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7. If both MppDL and Mmax exist and Mpppp < Mmax then all intermediate valnes of 
M bonnded by the cnrves M = Mpppp and M = Mmax wonld give positive potential 
snpersolitons. In fact, for < M < Mpppp, one can get a seqnence of PPSWs 
converging to the positive potential donble layer solntion at M = Mpppp whereas for 
MppDL < M < Mmax, One can get positive potential snpersolitons. For Mpppp < 
M < Mmax, amplitnde of the snpersoliton increases with increasing valnes of M bnt for 
the present system the amplitndes of the snpersolitons are bonnded and the maximnm 
amplitnde of the snpersoliton can be obtained at M = Mmax- As a resnlt, we get two 
types of PPSWs separated by a PPDL and the existence of hrst type PPSWs are 
restricted by Me < M < Mpppp whereas the existence of second type PPSWs (i.e., 
positive potential snpersolitons) are restricted by Mpppi < M < Mmax- A hnite jnmp 
between the amplitndes of two types of PPSWs separated by a PPDL can be easily 
verihed by plotting V (0) against 0 or by plotting 0 against 

8. We have nsed the following notations to interpret the different solntion spaces : C - 
Region of coexistence of both NPSWs and PPSWs, N - Region of existence of NPSWs, 
P - Region of existence of PPSWs and S - Region of existence of Snpersolitons. 

Now we consider the compositional parameter spaces or solntion spaces with respect to the 
parameter fi for increasing valnes of p starting from p = 0, i.e., there is no positron in the 
system. 

Solution Space w.r.t p when p = 0 : Fignre[T] shows the compositional parameter space 
or solntion space with respect to the parameter p for p = 0, i.e., there is no positron in the 
system. From this hgure, we have the following observations: 

• The system snpports NPSW for all M > Mc- 

• The system does not snpport any NPDL. 

• The system does not snpport any PPDL. 

• According to the dehnition of snpersolitons, the system does not snpport snpersoliton 
of any polarity becanse the system does not snpport donble layer solntion of any 
polarity. 
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• PPSWs start to exist if fi exceeds a critical value /ic and for all M lying within 
Me < M < Mmax. 

• The system supports coexistence of both NPSWs and PPSWs if p > Pc and Me < 
M < Mmax- For cTje = CTpe = 0.9, the value of pc is 0.145. 

Solution Space w.r.t p for 0 < p < 0.0001: Figure [2] shows the compositional parameter 
space or solution space with respect to the parameter p for p = 0.00001, i.e., a very small 
amount of positron is injected in the system. Here, the main important and interesting 
differences from the solution space with p = 0 are as follows: 

• The system supports NPDL along the curve M = M^pdl- 

• NPSWs are restricted for Me < M < Mnpdl, i-e., the existence region of NPSWs is 
bounded by the curves M = Me and M = Mnpdl and consequently, the system does 
not support any negative potential supersoliton. 

• The region of coexistence of NPSWs and PPSWs is bounded by the curves M = Me, 
M = Mmax and M = Mmpdl- 

• The system supports coexistence of NPDL and PPSW along the curve M = M^ppp 
bounded by the curves M = Me and M = Mmax- 

• As the system does not support any PPDL solution, the system does not support any 
positive potential supersoliton. 

For increasing values of p lying within 0 < p < 0.0001, the qualitative behaviour of the 
solution spaces remains unchanged, the only exception is that the region of existence of 
NPSWs decreases with increasing values of p for 0 < p < 0.0001. 

Solution Space w.r.t p for 0.0001 < p < 0.025 : If we further increase the concentration 
of positron from p = 0.0001 to p < 0.025, the solution space is not qualitatively same as 
solution space for 0 < p < 0.0001. The main difference in this solution space is the existence 
of PPDL solution along the curve M = Mpppp for a very small interval of p which has 
been shown in hgure [3] for p = 0.01. Again, it has also been observed that as p increases 
within 0.0001 < p < 0.025, the interval of p for the existence of PPDL solution increases. 
Again, the existence of PPDL solution implies the existence of PPSWs bounded by the 
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curves M = Me and M = Mppdl- From the solution space as given in figure [3l we can 
again conclude that the region of existence of NPSWs decreases with increasing values of 
p for 0.0001 < p < 0.025. This fact can be confirmed with the help of the solution space 
for p = 0.02 where we have seen a huge reduction in the existence region of NPSWs. The 
solution space with respect to the parameter p is given in figure IHfor p = 0.02. It is observed 
that the existence region of NPSWs ultimately collapses for p = 0.025. In fact, the existence 
region of NPSWs decreases for increasing values of positron concentration and finally, there 
exists a critical value of p for which the system does not support any negative potential 
solitary structure (NPSW or NPDL). This fact can be efficiently described by the figure 13 
The figure [5] is the solution space with respect to p for different values of p showing the 
existence regions of NPSWs bounded by the curves M = Me and M = Mpipph- But in case 
of the the existence region of PPSWs bounded by the curves M = Me and M = Mpppp, 
we have seen that the existence region of PPDLs increases for increasing values of positron 
concentration p ii p < p^^\ where p^^^ is a certain critical value of p whereas for p > p^‘^\ the 
existence region of PPDLs decreases for increasing values of positron concentration p and 
finally, there exists a critical value of p for which the system does not support any positive 
potential double layer. This fact can be efficiently described by the figure [6l The figure |6] is 
the solution space with respect to p for different values of p showing the existence regions 
of PPDLs. On the other hand, it can be easily verified that the existence region of PPSWs 
bounded by the curves M = Me and M = Mmax increases for increasing values of positron 
concentration p. For Uje = (Jpe = 0.9, the value of is 0.0215. Again from the solution 
spaces as given in figure [3] and figure 01 it is evident that the system does not support any 
positive or negative potential supersoliton. 

Solution Space w.r.t p for 0.025 < p < 0.07 : If we further increase the concentration 
of positron from p = 0.025 to p < 0.07, then the solution space is not qualitatively same as 
solution space for 0.0001 < p < 0.025. For p = 0.03, the solution space with respect to the 
parameter p is given in figure [71 Here we have the following observations: 

• The system does not support any NPSWs. 

• The existence region of PPSWs has increased. 

• There exist two cutoff values pa and p^ such that for pa < p < Pb, 'Vfe have Mpp^i < 
Mmax, be., in this interval of p, there exist PPSWs after the formation of double layers 
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if the mach number M is restricted by Mppdl < M < M^ax- In fact, we have two 
types of PPSWs: First Type - bounded by the curves M = Me and M = Mpppp and 
Second Type - bounded by the curves M = Mpppi and M = M^ax- Therefore, if 
these two types of PPSWs are separated by a PPDL, i.e., if the existence region or a 
portion of existence region of the total population of PPSWs (i.e., union of First and 
Second Type PPSWs) is separated by the curve M = Mpppp or a portion of the curve 
M = Mpppp, then for this particular portion of the curve M = Mpppp, we have a 
positive potential solitary wave after the formation of positive potential double layer. 
According to the property as mentioned by Das et al. [^, we have a jump type of 
discontinuity between the amplitudes of these two types of solitary waves just before 
and just after the formation of double layer and consequently, formation of positive 
potential supersoliton is conhrmed. Now, for the present problem, in the existence 
region of supersolitons, we have Mpppi < M < M^ax as supersoliton starts to exist 
just after the formation of PPDL and consequently, we have the bounded positive 
potential supersolitons. So, the region of existence of supersolitons is bounded by the 
curves fi = pa, fi = pb, M = Mpppi and M = Mmax- It is not difficult to hnd the 
values of /Xq and /if,. For Uje = CTp^ = 0.9, p = 0.03, the values of pa and /if, are 0.069 
and 0.081 respectively. 

• In Figure [8], the existence of positive potential supersolitons has been verihed by 

plotting V{(j)) against 0 and also by plotting 0 against ^ for a hxed /i lying within 
Pa < < l-^b and M = Mpppp + 0.0001 along with the other values of the parameters 

as shown in the hgure. In fact, we can take any p lying within pa < p < pb and any 
M lying within Mpppp < M < Mmax to verify the the existence of positive potential 
supersolitons. 

• It is observed that the existence region of supersolitons ultimately collapses. In fact, 
the existence region of supersolitons decreases for increasing values of positron con¬ 
centration and hnally, there exists a critical value of p for which the system does not 
support any positive potential supersoliton. This fact can be efficiently described by 
the hgure [9l 

From these observations, we can conclude that there must exist a critical value of p where 

the NPSW collapses and there must exist a critical value of p for which there exist two cutoff 
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values Ha and Hb such that for /!„ < /i < /x;,, we have Mppm < Mmax- 
Solution Space w.r.t /x for 0.07 < p < 0.117 : If we further increase the concentration 
of positron from p = 0.07 to p < 0.117, then the supersolitons disappear from the system 
and the existence region of PPDL is still decreasing. In this case, the system supports only 
PPSWs for all /x > 0. But there exists a cut off value pd of p such that the PPSWs are 
bounded by the curves M = Me and M = Mpppp for 0 < p < pd whereas for pd < p < Pt, 
the PPSWs are bounded by the curves M = Me and M = Mmax-, where pp is the physically 
admissible upper bound of p. Here the solution space with respect to the parameter p for 
p = 0.07 is given in hgure [TUI For the given values of cTje = 0.9 and ape = 0.9, the value of 
Pd is 0.018. 

Solution Space w.r.t p for p > 0.117 : If we further increase the concentration of positron 
starting from p = 0.117, then the PPDL disappears from the system but the system supports 
only PPSWs for all /x > 0 and these PPSWs are bounded by the curves M = Me and 
M = Mmax- For p = 0.2, the solution space with respect to the parameter p is given in 
hgure [m This solution space is not qualitatively same as solution space for p = 0.07. Here 
we see that the PPDLs disappear from the system but PPSWs exist for the entire range of 
p and the existence region of PPSWs is bounded by the curves M = Me and M = Mmax- If 
we further increase the concentration of positron then there is no qualitative change in the 
solution space, i.e., we get the similar type of solution space as obtained in case of p = 0.2 
but in any case the PPSWs are restricted by Me < M < Mmax and Mmax increases with 
increasing values of p, i.e., the region of existence of PPSWs restricted by Me < M < Mmax 
increases with increasing values of p starting from p = 0.117 and for the entire range of p. 

Next we shall consider the solution spaces with respect to p for some hxed value of p but 
for different values of ape- 

• For p = 0.01 and ape = 0.9, we have seen in hgure [3] that the system does not 
support any supersoliton but if the value of ape decreases then we have seen that 
NPDL disappears from the system and the positive potential supersoliton comes into 
the picture, i.e., there exists a critical value ap^J of ape such that the system does not 
support any supersoliton if ape exceeds its critical value apj and there exists a critical 
value ape^ of ape such that NPDL solutions come into the picture if ape > ape^ where 
cTpe^ > (Tpe^. For (Tje = 0.9, p = 0.01, the values of expe^ and ap^J are 0.5445 and 0.499 
respectively. 
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Figure [12] shows the solution space with respect to /i for p = 0.01 and ape = 0.4. This 
solution space shows the existence of positive potential supersolitons but nonexistence 
of NPDLs. In fact, for the decreasing values of ape starting from ape = 0.4, the system 
will always support supersolitons but such supersoliton disappears if ape exceeds the 
critical value a^}. 

Next, hgure [13] shows the solution space with respect to /r for p = 0.01 and ape = 0.5 
and we have seen that positive potential supersoliton disappears from the system. In 
hgure [TU the region of existence of positive potential supersolitons has been shown 
for different values of ape- From this hgure, we see that the region of existence of 
supersolitons increases with increasing values of ape- In fact, it has been observed that 
there exists a critical value ape of ape such that at ape = ap} the region of existence of 
supersolitons is maximum and after ape = apj the region of existence of supersolitons 
decreases very rapidly with increasing values of ape and hnally it disappears from the 
system and after the disappearance of the positive potential supersolitons, NPDLs 
come into the picture as we have seen earlier for ape = 0.9. For ate = 0.9, p = 0.01, 
the values of apj is very close to a^J. 

Other interesting results obtained from the solution spaces are the followings: 

— Region of existence of PPSWs bounded by the curves M = Me and M = Mmax 
decreases with increasing values of ape as shown in the hgure [15] 

— Region of existence of PPSWs bounded by the curves M = Me and M = Mpp^i 
increases with increasing values of ape as shown in the hgure [16] 

— It can also be checked that the region of existence of NPSWs bounded by the 
curves M = Me and M = Mpfppp increases with increasing values of ape- In fact, 
here we have qualitatively same type of solution spaces as given in hgure [5] The 
only exception is that the region of existence of NPSWs bounded by the curves 
M = Me and M = M^ppp increases with increasing values of ape instead of the 
fact that the region of existence of NPSWs bounded by the curves M = Me and 
M = Mjsipdl decreases with increasing values of p. 

— For small value of p and for very large value of ape {ape = 100, say), we have 
exactly same type of solution space as given in hgure [T] i.e., the solution space 


for p=0. This fact is consistent with the Energy integral (|2^ becanse it is simple 
to check that 


lim apeWp = —, 

(Tpe—>-00 <h 


and conseqnently the eqnation fl2^ can be written as 
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1 — p 1 — p 

Now, for small valnes of p, we can set 1 — p 1 and the above expression of W (</>) 
assnmes the following form 
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The eqnation fl6TD is free from the parameters associated with the positrons and 
conseqnently, for this case, there is no effect of positrons on the solitary strnctnres 
of the DIA waves. On the other hand, if p+p < 1, then we have the qnalitatively 
similar type of solntion space as given in hgnre [H 


From the solntion space with respect to p for different valnes of ape (actnally, for increasing 
valnes of ape), we have seen that the NPDL comes into the pictnre after the disappearance 
of positive potential snpersolitons whereas from the solntion space with respect to p for 
different valnes of p (actnally, for increasing valnes of p), we have seen that the positive 
potential snpersoliton comes into the pictnre after the disappearance of NPDL. Therefore, 
for the solntion space with respect to p, we have seen that it is impossible to get simnltaneons 
occnrrence of positive potential snpersolitons and NPDLs for any arbitrarily hxed valnes of p 
and ape- This fact can be easily nnderstood from the solntion space with respect to p. Fignre 
[T7I shows the solntion space with respect to p for hxed valnes of p and ape as mentioned in 
the hgnre. Althongh this solntion space shows the simnltaneons occnrrence of NPDLs and 
positive potential snpersolitons bnt if we restrict the valne of p a.t p = pk (a hxed valne of 
p), then from the hgnre it is clear that the simnltaneons occnrrence of NPDL and positive 
potential snpersoliton is not possible. In other words, if we move the solntion space along 
the cnrve parallel to the cnrve M = Me then from the solntion space as given in hgnre [T71 we 
see that the positive potential snpersoliton comes into the pictnre after the disappearance 
of all NPDLs. Again, we have seen the coexistence of NPSWs and PPSWs for an interval 
of p bonnded by the cnrves M = Me, M = and M = M^ax only when the positive 
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potential supersoliton disappears from the system. But this case occurs only when /i exceeds 
a critical value. This qualitative behaviour has been shown in the £gure[T8l Figure [18] shows 
the solution space with respect to p for fixed values of p and ape as mentioned in the hgure. 
From this solution space we see that positive potential supersoliton disappears from the 
system and at the same time coexistence of both positive and negative potential solitary 
waves take place for a dehnite interval of p. If we further increase the value of p then we 
see that NPDLs disappear from the solution space and for the entire interval of p, we have 
only PPSWs bounded by the curves M = Me and M = M^ax- This type of solution space 
is qualitatively same as the solution space with respect to p as shown in hgure [TTJ 


VI. SUMMARY & DISCUSSIONS 

In the present paper, a thorough investigation on the existence regions of the different 
dust ion acoustic solitary structures has been made in an unmagnetized dusty plasma con¬ 
sisting of negatively charged dust grains, adiabatic warm ions, and isothermally distributed 
electrons and positrons with the help of the compositional parameter space showing the 
nature of existence of different solitary structures, giving a special emphasis on the nature 
of existence of supersolitons. If there exist solitons after the formation of double layer of 
same polarity for increasing values of the Mach number then the solitons are known as su- 
persolitons and if the amplitude of any supersoliton is bounded then the supersolitons are 
known as bounded supersolitons. The Sagdeev potential technique has been used to study 
the arbitrary amplitude supersolitons. A general theory has been given on the existence of 
bounded supersolitons. If there exists a parameter regime for which the positive (negative) 
potential double layer is unable to restrict the occurrence of positive (negative) potential 
solitary waves then in this region of the parameter space, there exist positive (negative) po¬ 
tential solitary waves after the formation of positive (negative) potential double layer, and 
consequently, formation of positive (negative) potential supersoliton is confirmed. On the 
other hand, if the existence region or a portion of existence region of the total population 
of PPSWs (NPSWs) is separated by the curve M = Mpp^i (M = M^ppp) or a portion 
of the curve M = Mpppp (M = Mnpdl), then for this particular portion of the curve 
M = Mpppp {M = Mpfppp), we have a positive (negative) potential solitary wave after the 
formation of positive (negative) potential double layer and as a result, we have a jump type 
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discontinuity between the amplitudes of two types of solitary waves just before and just after 
the formation of double layer and consequently, formation of supersoliton is conhrmed. In 
the present problem, we have found positive potential bounded supersolitons. The present 
dusty plasma system does not support any negative potential supersoliton. Not only /i and 
p but also the parameter cXpe plays an important role for the formation of bounded positive 
potential supersolitons. Other than the formation of bounded positive potential supersoli¬ 
tons, we have found many interesting results regarding the existence of solitary structures 
from the qualitatively different compositional parameter spaces or solution spaces showing 
the nature of existence of different solitary structures with respect to any parameter of the 
system. One of the interesting observation is that the NPDLs occur in a deleted right neigh¬ 
bourhood of p = 0 (see the solution spaces given in the hgure [IT] and hgure HSj) whereas 
PPDLs exist in a deleted right neighbourhood of p = 0 (see the solution spaces given in the 
hgure |3l hgure H] hgure [7| and hgure [10]). These observations are very much consistent with 
the physical property for the formation of double layer which states that for the formation 
of NPDL (PPDL) solution at a point in the compositional parameter space the negative 
(positive) potential (absolute value) must dominate the positive (negative) potential in a 
neighborhood of that point of the compositional parameter space and the potential diher- 
ence must be maximum there. Another interesting observation is that for small value of 
positron concentration there is no qualitative change in solitary structures of DIA waves if 
very hot positrons are injected in the system. This fact is also consistent with the Energy 
Integral (l22ll . Now, we want to summarize the main observations from the qualitatively 
different compositional parameter spaces with respect to p for different values of p and ape 
and also with respect to p for different values of /i. 

1. When there is no positron in the system, the system supports NPSW solution whenever 
M > Me, i.e., the existence region of NPSWs is unbounded above whereas the existence 
region of PPSWs is bounded by the curves M = Me and M = Mmax- 

2. Whenever a small amount of positron is injected in the system, the system supports 
NPDL solution along the curve M = M^pol and consequently, the existence region 
of NPSWs has become bounded by the curves M = Me and M = Mjsipdl- The 
existence region of NPSWs decreases for increasing values of positron concentration 
and hnally, there exists a critical value of p for which the system does not support 
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any negative potential solitary structures (NPSWs and NPDLs). This fact can be 
efficiently described by the hgure [5l The hgure [5] is the solution space with respect 
to /i for different values of p showing the existence regions of NPSWs bounded by 
the curves M = Me and M = If the negative potential solitary structures 

leave the system for a particular value of p then it will never come back for further 
increment in positron concentration. 

3. When there is no positron in the system, the system does not support any PPDL 
solution. If we increase p from p = 0, then the system will support the PPDL solution 
whenever p exceeds a cut off value pi (say). Then the existence region of PPSWs 
bounded by the curves M = and M = Mpppp will increase for increasing values 
of p but here also there exists a cut off value p 2 (say) such that the existence region 
of PPSWs bounded by the curves M = Me and M = Mpppp increases for increasing 
p for Pi < p < p 2 whereas the existence region of PPSWs bounded by the curves 
M = Me and M = Mpppi will decreases for increasing p for P 2 < P < Ps and at 
p = P 3 , the PPDL solution disappears from the system. This fact can be efficiently 
described by the hgure El The hgure El is the solution spaces with respect to p for 
diherent values of p showing the existence regions of PPSWs bounded by the curves 
M = Me and M = Mpppp. If the positive potential double layer leaves the system 
for a particular value of p then it never comes back for further increment in positron 
concentration. But there always exist positive potential solitary waves bounded by 
the curves M = Me and M = Mmax- So, although the PPDL solution disappears from 
the system for p = pa, there always exist PPSWs bounded by the curves M = Me and 
M = Mmax- In fact, for p > ps, the system supports only PPSWs for any value of p. 

4. From the previous discussions we have seen that if p exceeds the cut oh value pi 
but the value of p is restricted by the inequality p < Ps then we have two types of 
PPSWs: Type - I : PPSWs bounded by the curves M = Me and M = Mpppi which 
is same as First type PPSWs as dehned earlier and Type - II : PPSWs bounded 
by the curves M = Me and M = Mmax which is not same as Second type PPSWs 
(see the solution space as given in hgure El hgure 01 hgure [71 hgure 001) • Now if 
these two types of PPSWs are separated by a PPDL, i.e., if the existence region or 
a portion of existence region of the total population of PPSWs (i.e., union of Type 
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- I and Type - II PPSWs) is separated by the curve M = Mppdl or the portion of 
the curve M = Mpppp, then for this particular portion of the curve M = Mpppp, 
we have a positive potential solitary wave after the formation of double layer (see the 
region marked as ‘S’ of the solution space as given in hgure [7]). Now according to 
the property mentioned by Das et al. [70|, we have a jump type discontinuity between 
the amplitudes of two types of solitary waves just before and just after the formation 
of double layer and consequently, formation of supersoliton is conhrmed (see hgure 
[8](d)). Now, the upper bound of PPSWs is either Mmax or Mpppp. Obviously, in the 
existence region of positive potential supersolitons, we have Mpppp < Mmax as positive 
potential supersoliton starts to exist just after PPDL. Therefore, for the bounded 
positive potential supersolitons, the following inequality holds good: Mpppi < M < 
Mmax, where M is the Mach number corresponding to a positive potential supersoliton, 
i.e., a PPSW just after the formation of positive potential double layer. Figure [9] 
shows that the existence region of the second type of PPSWs (i.e., positive potential 
supersolitons) decreases for increasing values of p and hnally disappears from the 
system if p exceeds a certain cut off value. 


5. For hxed ape, region of existence of PPSWs bounded by the curves M = Me and 
M = Mmax increases with increasing values of p. 

6. For hxed ape, region of existence of PPSWs bounded by the curves M = Me and 
M = MpppL increases for increasing p for pi < p < p 2 whereas it decreases for 
increasing p for p 2 < p < Ps- 


7. For hxed ape, region of existence of NPSWs bounded by the curves M = Me and 
M = M]\fpp)p decreases with increasing values of p. 


8. For hxed ape, region of existence of supersolitons increases with decreasing values of p 
up to a certain critical value of p and then it decreases and hnally it disappears from 
the system. 


9. For hxed p, region of existence of PPSWs bounded by the curves M = Me and M = 
Mmax decreases with increasing values of ape- 

10. For hxed p, region of existence of PPSWs bounded by the curves M = Me and M = 
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MppDL decreases with increasing values of ape. 

11. For fixed p, region of existence of NPSWs bounded by the curves M = Me and 
M = Mjqppp increases with increasing values of ape- 

12. For fixed p, region of existence of positive potential supersolitons increases with in¬ 
creasing values of ape up to a certain critical value of ape and then it decreases and 
hnally it disappears from the system. 

13. From the solution space with respect to p, we have seen that for a particular p, it is 
not possible to get coexistence of NPSWs or NPDL and positive potential supersoliton 
whereas from the solution space with respect to p, we have seen that for a particular 
p, it is not possible to get coexistence of NPSWs or NPDL and positive potential 
supersoliton. 

14. From hgure [IZl it is clear that for a particular value of p, occurrence of NPDL does 
not imply the occurrence of positive potential supersoliton and conversely, occurrence 
of positive potential supersoliton does not imply the occurrence of NPDL, i.e., in 
any case, simultaneous occurrence of NPDL and positive potential supersolton is not 
possible. In fact, we have seen earlier that positive potential supersoliton comes into 
the picture after the disappearance of NPDLs for the increasing ( decreasing ) values 
of p {ape) whereas NPDL comes into the picture after the disappearance of positive 
potential supersolitons for the increasing (decreasing) values of ape {p). 

15. NPDL exists for small values of p whereas PPDL exists for small values of p. 

16. For small value of p and for very large value of ape (ape = 100, say), we have exactly 
same (qualitatively similar) type of solution space as given in hgure [H i.e., the solution 
space for p=0. This fact is consistent with the Energy integral ([22]). 

To conclude, we want to mention that the investigations made in the present paper 
through the qualitatively different compositional parameter spaces or solution spaces are 
helpful for the understanding of different nonlinear electrostatic structures of the present 
plasma system for the physical parameters of our interest. Our investigations are useful 
in understanding the properties of nonlinear dynamics of ion acoustic perturbations that 
may appear in astrophysical plasmas, such as those in the early universe and active galactic 
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nuclei. In the present paper, we have not investigated the nonlinear behaviour of DIA waves 
when U = cd ^ W"{Q) = 0 ^"(0) = 0 M = Me, i.e., when the velocity of the wave 

frame is equal to the linearized velocity of the DIA wave for long wave length plane wave 
perturbation. The study of the dust ion acoustic solitary structures at the acoustic speed 
{cd) in such plasmas is under consideration and the results shall be reported elsewhere. 
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1.07 



FIG. 1. Compositional parameter space with respect to fi ioi p = 0, aie = 0.9 and ape = 0.9 



0 jj^ ^ 0.5 


FIG. 2. Compositional parameter space with respect to p for p = 0.00001, aie = 0.9 and ape = 0.9 


39 





40 







a. = 0 . 9,0 = 0.9 
le pe 



FIG. 5. Compositional parameter space with respect to n showing the negative potential solitary 
structures for different values of p 

a. = 0.9, a = 0.9 
le pe 



FIG. 6. Compositional parameter space with respect to p showing the positive potential double 
layers and positive potential solitary waves restricted by positive potential double layers for different 
values of p 
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0 0.15 


FIG. 7. Compositional parameter space with respect to fi for p = 0.03, crje = 0.9 and ape = 0.9 





FIG. 8. V{(p) is plotted against (/> for a value of p lying in the interval //„ < M < /r;, for p = 0.03: 
(a) M = MppDL — 0.0001 (b) M = Mpppp (c) M = Mpppp + 0.0001. In (d), cj) is plotted against 
^ for the same values of p and p for M = Mpppp — 0.0001 (dash-dot) and M = Mpppp + 0.0001 
(solid). 
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FIG. 9. Compositional parameter space with respect to ^ showing the existence region of super- 
solitons for different values of p 



0 0.018 0.1 


FIG. 10. Compositional parameter space with respect to p for p = 0.07, aie = 0.9 and ape = 0.9 







FIG. 11. Compositional parameter space with respect to lor p = 0.2, aie = 0.9 and ape = 0.9 



0 0.2 


FIG. 12. Compositional parameter space with respect to p for p = 0.01, aie = 0.9 and ape = 0.4 
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FIG. 13. Compositional parameter space with respect to n for p = 0.01, aie = 0.9 and ape = 0.5 
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FIG. 14. Compositional parameter space with respect to /x showing the existence region of super- 
solitons for different values of ape 
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0 0.2 


FIG. 15. Compositional parameter space with respect to ^ showing the positive potential solitary 
waves bounded by M = Mmax for different values of ape 



FIG. 16. Compositional parameter space with respect to n showing the positive potential double 
layers and positive potential solitary waves restricted by positive potential double layers for different 
values of ape 
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FIG. 17. Compositional parameter space with respect to p for ^ = 0.05, aie = 0.9 and ape = 0.9 
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FIG. 18. Compositional parameter space with respect to p for p = 0.15, aie = 0.9 and ape = 0.9 
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